Anomalous Levy decoherence 
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PACS numbers: 05.40.Fb, 03.65.Yz, 05.40.-a 

First a subject of mathematical investigation ([10, 
Levy stable processes have nowadays fully entered the 
realm of physics y,Q] . It is well known that usual Brow- 
nian motion is a Gaussian stochastic process [g|. This 
is a consequence of the Gaussian law of large numbers, 
or central-limit theorem. In the 1930s, the central-limit 
theorem has been extended by Levy and others to ran- 
dom variables with infinite second moment. More pre- 
cisely, the probability distribution of these random vari- 
ables has the form, 
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We investigate the decoherence of a small quantum system weakly coupled to a complex, chaotic 
environment when the dynamics is not Gaussian but Levy anomalous. By studying the time de- 
pendence of the linear entropy and the damping of the interference of two Gaussian wave packets in 
the Wigner representation, we show that the decoherence time for a quantum Levy stable process 
is always smaller than for Gaussian diffusion. 



Decoherence manifests itself in the dynamical destruc- 
tion of quantum interferences as a result of the interac- 
tion with the environment [pL4|| . This progressive transi- 
tion from quantum to classical has been experimentally 
observed in high Q microwave cavities |jj| ■ Theoretically, 
decoherence has been intensively investigated in the case 
of normal (Gaussian) diffusion [[l6|-[l8| but not so far for 
anomalous Levy diffusion. 

A simple way to evaluate the decoherence time r^ is 
to compute the linear entropy S a (t) = 1 — tr p 2 (t), which 
gives a measure of the purity of the system |19| , |20|| : The 
linear entropy is equal to zero for a pure state and in- 
creases to one when the state of the system evolves to a 
mixture. Here p(t) is the reduced density operator of the 
system. Alternatively, one can consider the superposition 
of two Gaussian wave packets in the Wigner representa- 
tion and look at the damping of the interference fringes 
pl[ . In the following we study the time dependence of 
both the linear entropy and of the interference term of 
the Wigner function. We put a special emphasis on the 
value a = 1 where analytical results can be obtained. 
What we find is that Levy decoherence is always faster 
than Gaussian decoherence for any value of the stability 
index a. 

Starting point of our analysis is a master equation 
for a system coupled to a chaotic environment derived 
by Kusnezov et al. from a microscopic random-matrix 
model ||. In this approach the coupling to the non- 
homogeneous background is characterized by a spread- 
ing width F^ and a spatial correlation function G(x/xq) 
with an environment characteristic correlation length xq 
(see Ref. ||] for details). In the limit of high (actually 
infinite) temperature, the quantum master equation for 
a free particle of mass M is given by 
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0<a<2. (1) 



A (symmetric) Levy stable distribution as given by ([fl) 
can be regarded as an extension of the Gaussian distribu- 
tion to which it reduces when a = 2. It should be stressed, 
however, that the Gaussian distribution is the only stable 
law having a finite variance. A Levy process may then 
be defined as a non-Gaussian generalization of Brownian 
motion obeying Levy stable statistics. The hallmark of 
this form of anomalous diffusion is the presence of long 
jumps (the so-called "Levy flights") which are due to 
the asymptotic power-law tail of a stable distribution, 
£„(x) ~ l/|x| Q+1 . As a consequence, the mean-square 
displacement and the mean kinetic energy of a Levy par- 
ticle are divergent |6|. Levy dynamics often results from 
the interaction with a complex, non-homogeneous envi- 
ronment. Examples are porous or disordered media |7]] 
or chaotic heat baths M. 

The first experimental observation of a Levy stable 
process has been reported ten years ago in micelle sys- 
tems M. Since then, Levy-like diffusion has for in- 
stance been seen in two-dimensional rotating flows (HI , 
in porous glasses IO] and, more recently, in ion lattices 
Jr2[ and in subrecoil laser cooling [Of. The ion lattice 
experiment is of particular interest since it reports a di- 
rect measurement of the divergence of the mean kinetic 
energy of the particle. Moreover, it is worthwhile to no- 
tice that experiments have now started to study anoma- 
lous Levy dynamics of microscopic systems (atoms, ions) 
where quantum effects are likely to play a role. This mo- 
tivated us to examine the decoherence of quantum Levy 
stable motion. 
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When the environment correlation length xq is much 
larger than the separation Ax = \x — y\, the correlator 
in Eq. (pi) can be expanded as p3] 
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The case a = 2 corresponds to Gaussian diffusion, while 
for a faster spatial decorrelation, <a< 2, the diffusion 
is Levy anomalous [g 23 j. The (generalized) diffusion co- 
efficient is further defined as D a — T^/fixy). It is also 
worthwhile to emphasize that it is the last term in the 
master equation (||), which results from the accumula- 
tion of phase fluctuations of the environment, which is 
mainly responsible for the decoherence effect |Q. We 
already know that the alteration of this noise term when 
passing from Gauss to Levy statistics has a dramatic ef- 
fect on the diffusion properties of the system, it is thus 
of high interest to examine the effect of the modification 
of this term on decoherence. 

We begin by solving the master equation (||) for p(t). 
For this purpose it is convenient to use the (k, A) repre- 
sentation of the density matrix p4| , 



p(k,A,t) = 



dqe tkq p{q+- 1 q-- 1 t) , 



(4) 



instead of the coordinate representation p(x, y, t). Here, 
we have introduced the center of mass and relative coor- 
dinates, q = (x + y)/2 and A = x — y. In this represen- 
tation, the master equation reads 

|p(*,A,t) = {A^_2) a |Ar}p(*,A,t). (5) 

This first order partial differential equation can easily be 
solved, for instance by the method of characteristics . We 
find 
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where p(k, A, 0) is the reduced density operator at t = 0. 
For an initial Gaussian pure state, p(x,y 1 0) — (p(x)if*(y) 
with ip(x) = exp[— x 2 /46 2 ]/(27t6 2 ) 1//4 , it has the form 



p(k,A 7 0) — exp 
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It is further interesting to consider the probability distri- 
butions for position and momentum that correspond to 
the solution (|6|). In the limit of long times, we have 
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dk e- lkx p{k, A = 0, t) ~ £«« (x) , (8) 
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with u(t) = D a t 1+a /(l + a)M a and v(t) = D a t. The 
result (ph, which holds for 1 < a < 2, has already been 



obtained in Ref. ||. Notice that both the position and 
the momentum distributions of the system asymptoti- 
cally approach a Levy stable law (in the short time limit 
they are given by Gaussians). In particular, this means 
that for a — 1 , p(p, p, t) has the form of a Lorentzian. It 
is worth mentioning that a recent experiment confirmed 
that the momentum distribution of atoms cooled by se- 
lective coherent population trapping is very close to a 
Lorentzian distribution p5fl . 

Linear entropy. Let us now calculate the linear en- 
tropy S a (t) of the system. We use the expression (|6|) 
for p(k,A,t) together with the initial condition (R). For 
short times, we then obtain 
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Consequently, the system gets entangled with the en- 
vironment within a decoherence time of the order of 
t d ~ [(rV^XV^o)"]" 1 , where b is the width of the ini- 
tial wave packet. 

In order to obtain more tractable expressions, we 
will work in the following in the limit h — ► 0, where 
D a ((Ax) a ) is large |0|. In this "semiclassical" limit, 
the time evolution of the density operator is dominated 
by the last term in Eq. (||). It is then straighforward to 
write down an approximate solution of the master equa- 
tion in the form, 



p{x,y,t) = p(x,y,0) exp 
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We observe that the off-diagonal matrix elements decay 
exponentially on a time scale r^ = [(r^/?i)(Aa;/a;o) Q! ] _1 . 
Note that this expression of the decoherence time is in 
agreement with the previous result obtained from the 
short time behavior of the linear entropy. Let us briefly 
discuss the validity of this approximation. To this end, 
it is instructive to write the approximate solution (|ll]) of 
the master equation in the (k, A) representation which 
yields 



p(k,A,t) = p(k,A,0)exp[-D a \A\ a t} 



(12) 



This expression is easily seen to be a short time approx- 
imation of the exact solution (0) [see also Fig. (Eh]. 
Superposition of wave packets. We now turn to the 
study of the effect of decoherence on the interference of 
two Gaussian wave packets in the Wigner representation. 
The Wigner transform of the density matrix is defined as 
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We mention that the (pseudo) phase space distribution 
W(q,p,t) is also related to the (k, A) representation by 
a (double) Fourier transform |l(|. Applying the Wigner 
transform to the master equation (0) (retaining again 
only the last term) we arrive at a fractional diffusion 
equation in momentum, 



lw(q,p,t)=D a ^W(q,p,t) 



(14) 



Here the Riesz fractional derivative is defined through its 
Fourier transform as EGlET! , 
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The solution of the fractional equation (14) can for in- 
stance be found in J6J. Next we consider an initial state 
of our system that consists of the linear superposition of 
two Gaussian wave packets separated by a distance 2a: 
tp(x) = ciexp[-(a; - a) 2 /4b 2 ] + c 2 exp[-(a; + a) 2 /4& 2 ]. 
The Wigner function corresponding to this initial state 
can be written as a sum of two classical contributions plus 
a quantum (interference) term: W — W\ + Wi + Wi nt . 
For the value a= 1, analytical evaluation of the Wigner 
function W(q,p, t) is possible. This leads to 
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^{c\c 2 f{a) + c l c* 2 f{-a)} . 
The function /(a) is given by the expression 



(16) 
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(18) 



Erf(x) = 2/y/ir f Q exp[— t 2 ]dx being the Error function. 
Moreover, in the limit where the two wave packets are 
well separated (a > d), the time dependence of the in- 
terference term can be simplified to become 



W int (q, 0,£) ~ a/ — exp 
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Hence the interference between the two superposed Gaus- 
sian wave packets is exponentially destroyed by decoher- 
ence on a time r^ = [(r^/?i)2a/a;o] _1 , where 2a is the 
distance between the two wave packets. 
Discussion. Let us now discuss the results that we have 
obtained so far. We have found that for a quantum Levy 



process with exponent a, the suppression of coherence 
over a distance Ax takes place exponentially on a time 
scale rg ~ [(T^^^Ax/xq)"}' 1 , Ax being either the 
width of the single wave packet in the computation of 
the linear entropy or the separation of the two packets in 
the evaluation of the Wigner function. The decoherence 
time t^ is hence inversely proportional to the spreading 
width r^ and also inversely proportional to the ath power 
of the ratio of the distance Ax over the environment cor- 
relation length xq . When compared to the Gaussian case 
a = 2, we find that t^/t^ ~ (Ax/xo) 2 ~ a , independent 
of the factor T^ /%. 
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FIG. 1. Linear entropy S a (t) as a function of time for dif- 
ferent values of the stability index a: The curves (a) to (d) 
are obtained with the exact solution (H) of the master equa- 
tion for increasing values of a, a = 0.5, 1, 1.5, 2. The curve 
(e) corresponds to the approximate solution (\12\j in the case 



In the present analysis the environment correlation 
length xq is much larger than any typical system sep- 
aration Ax, so that we have Ax/xq <C 1. Since for Levy 
stable motion a < 2, we are therefore lead to the conclu- 
sion that Levy decoherence is always (for any exponent 
a) faster than Gaussian decoherence (or in other words 
its decoherence time is always smaller) . The dependence 
of the decoherence time on the value of the stability in- 
dex is shown in Fig. (|lj), where the linear entropy S a (t) 
is plotted as a function of time for different values of the 
parameter a. We see that the smaller the value of a, the 
faster the decoherence for a Levy process. For instance, 
for a ratio Ax/xo = 10~ 3 , Levy decoherence is three 
orders of magnitude faster than Gaussian decoherence 
when a — I. Levy anomalous diffusion hence appears, 
in this sense, to be "more classical" than its Gaussian 
counterpart. 

This result can be interpreted as follows. In the 
Langevin picture of Brownian motion, the coupling to the 
environment manifests itself in the addition of two forces 
to Newton's equation: first, a friction force which dissi- 
pates energy away from the system and, second, a fluctu- 
ating stochastic force which randomly injects energy back 
into the system. For Levy stable motion, the variance of 



this stochastic force is divergent which means that the 
environment randomly supplies an infinite amount of en- 
ergy to the system (this is, by the way, the physical origin 
of the Levy flights) pq] . Thus, in this case, the coupling 
to the surroundings has a much stronger influence on the 
evolution of the system and, in turn, decoherence, like 
diffusion, is greatly enhanced. Note, however, that in 
contrast to diffusion, Levy decoherence does not involve 
divergent quantities. 

In summary, we have investigated the decoherence 
properties of a small quantum system weakly coupled to 
a chaotic background when its dynamics exhibits Levy 
type anomalous behavior. We have determined the de- 
coherence time, first, by computing the linear entropy of 
the system and, then, by examining the decay of interfer- 
ence fringes of two superposed Gaussian wave packets in 
the Wigner representation. In comparison to the Gaus- 
sian case, we have found that Levy decoherence is always 
faster. More precisely, the deviation of the decoherence 
time from its Gaussian value is larger the smaller the 
value of the characteristic exponent a. 
This work was supported by the RTD network Nanoscale 
Dynamics, Coherence and Computation. We thank D. 
Kusnezov and C. Lewenkopf for discussion. 
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